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Gaussianity of temperature fluctuations in the Cosmic Microwave Background(CMB) implies that 
the statistical properties of the temperature field can be completely characterized by its two point 
correlation function. The two point correlation function can be expanded in full generality in the 
bipolar spherical harmonic(BipoSH) basis. Looking for significant deviations from zero for Bipolar 
Spherical Harmonic(BipoSH) Coefficients derived from observed CMB maps forms the basis of the 
strategy used to detect isotropy violation. In order to quantify "significant deviation" we need to 
understand the distributions of these coefficients. We analytically evaluate the moments and the 
distribution of the coefficients of expansion(^4^^), using characteristic function approach. We show 
that for BipoSH coefficients with M = an analytical form for the moments up to any arbitrary 
order can be derived. For the remaining BipoSH coefficients with M 0, the moments derived using 
the characteristic function approach need to be supplemented with a correction term. The correction 
term is found to be important particularly at low multipoles. We provide a general prescription for 
calculating these corrections, however we restrict the explicit calculations only up to kurtosis. We 
confirm our results with measurements of BipoSH coefficients on numerically simulated statistically 
isotropic CMB maps. 
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I. INTRODUCTION 



Cosmological model building has been usually pursued under the assumption that the universe is homogeneous 
and isotropic. Statistical isotropy of CMB implies statistical expectation values of the temperature fluctuations are 
preserved under rotations of the sky. The CMB data is one of the cleanest observations and it is only reasonable to 
search for weak violations of statistical isotropy in the CMB maps. The source of this deviation remains mysterious 
so far. However, with present as well as future maps providing more and more information it is important to study 
deviations from standard statistics. 

Random temperature fluctuations in the CMB arc believed to be Gaussian, these fluctuations can be completely 
characterized by specifying the two point correlation function. 

The two point correlation function is most generally expanded in the bipolar spherical harmonic(BipoSH) basis 
to test the violations of isotropy in the CMB temperature and polarization maps. This formalism was developed by 
Hajian and Souradeep [H-@ and is such that for an isotropic sky all BipoSH coefficients, Af^ except A^j 3 vanish on 
an average. These expansion coefficients have been used to parametrize several kind of statistical isotropy violations 
0413 an d was adopted by the WMAP team 11] to search for violations of isotropy in the WMAP data. Although, 
these coefficients were primarily introduced to study statistical isotropy violation, they have found various other 
applications [l^ - [l^ |. 

Specifying all the moments of a distribution completely characterize the distribution. In this paper, we derive 
analytical expressions for the moments of the distribution of the BipoSH coefficients using the characteristic function 
approach. BipoSH coefficients are linear combinations of elements of the harmonic space covariance matrix. The 
independence of the terms in the linear combination for the BipoSH coefficients with M = ensures that the 
characteristic function encodes complete statistical information. For the remaining BipoSH coefficients with M =/= 0, 
the characteristic function method partially works due to the presence of non linear correlations among terms in the 
linear combination. To account for these non linear correlations we supply a correction term to the moments derived 
using the characteristic function method. We test these analytical results against simulations. We simulate the CMB 
maps using the widely used HealPix 15] package. 

This paper is organized as follows. In section [TT] we briefly discuss the BipoSH formalism introduced by Hajian 
and Souradeep. In section Hill we discuss the characteristic function approach which is extensively used to derive the 
moments of the distribution of the BipoSH coefficients. In section HVl we present the analytical expressions derived 
for the various BipoSH coefficients. The details of these calculations and a detailed discussion on the characteristic 
function approach can be found in the appendices. We conclude with a discussion of our results in section IVl 



II. BIPOSH FORMALISM 



The isotropic CMB sky is fully characterized by specifying the four angular power spectra CT , Cf , Cf and C\ T . 
These are the Legendre polynomial coefficients of the corresponding two point correlation function defined in the 
following manner, 

C xx (h h h 2 ) = C xx {h x .h 2 ) = V ^±V*P^i • n 2 ) . (1) 
In what follows we drop the 'XX' label for notational brevity. 

If the CMB sky is not assumed to be isotropic then two point correlation function in general will depend on the 
directions hi and h 2 . Hence, the bipolar spherical harmonic basis form a very natural basis in which the CMB two 
point correlation function can be expanded, 

C(fti,n2)= Y, A il{ Y h(ni)®Y l2 (h 2 )} L M, (2) 

h,h,L,M 

where Af¥ are BipoSH coefficients and {Yi 1 {ni)^Yi 2 (h 2 )}LM are bipolar spherical harmonics [l6j. BipoSH functions 
are irreducible tensor product of two spherical harmonics with different arguments, they form an orthonormal basis 
on S 2 x S 2 for different sets of l±,l 2 , L, M . Their transformation properties under rotations are similar to spherical 
harmonics and can be expressed as, 

{Y lx {h 1 ) ® Y h (h 2 )} LM = ]T c^^Aihx) Y% 2 (h 2 ), (3) 

minis 

where C^^fiJama are Clebsch- Gordon coefficients. These indices satisfy triangularity conditions \li — l 2 \ < L < l x + 1 2 
and mi + m 2 = M. 
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The BipoSH coefficients can be shown to be iinear combinations of off-diagonal elements of the harmonic 
space covariance matrix |l| , 

ra\ra 2 

where a; m 's are the spherical harmonic coefficients of the CMB maps. An unbiased estimator of BipoSH coefficients 
can be defined in terms of the spherical harmonic coefficients of the CMB maps, 

A hh = a hm 1 a,l 2 m 2 Gl J i 1 m, 1 l2m2 • (^) 

mim2 

It can be proved that for an isotropic CMB sky the expectation value of all the BipoSH coefficients, except the 
isotropic angular power spectrum Aff = (— l) ; C/V2i + 1 vanish [H Hoj. 



III. CHARACTERISTIC FUNCTION METHOD 



We investigate the statistical properties of the real and imaginary parts of complex coefficients obtained in the 
BipoSH representation of the CMB two point correlation function. To arrive at the moments of BipoSH coefficients, 
which are linear combinations of covariance matrix elements (see eq |C3p , we adopt the characteristic function approach 
which is widely used in statistics 17| . 

The characteristic function of any random variable completely defines its probability distribution fl8l | . It is defined 
in the following manner, 

<px{t) = E[e itx \ teR. (6) 

Consider a random variable defined in the following manner, 

n 

Z n = y^ j aiXj ) (7) 
»=i 

where Oj's are constants and X^s are independent random variables which are not necessarily identically distributed. 
By independence we imply that, all the higher order correlations between the terms appearing in the linear 
combination vanish, 



{X n . X *m )=Q V 7l : TYL (i^j). (8) 

The characteristic function method is particularly useful in arriving at the statistics of such random variables. The 
characteristic function of Z n will just be the product of the characteristic function of the individual terms contributing 
to the linear sum, 

(*) = <PXi (ait)<fix 2 (a 2 t) ip Xn (a n t) . (9) 

If the terms involved in the linear combination are not independent then the characteristic function will not take 
up the simple form given above. The cumulant generating function is defined as the logarithm of the characteristic 
function, 

g z (t)=\og[<p z (t)}. (10) 

The cumulants can be obtained by taking derivatives of the cumulant generating function and evaluating them at 
zero 

K n = i n g n z (t)\t=o- (11) 

Given the cumulants, it is straightforward to arrive at the moments of the distribution. The explicit relationships 
between cumulants and central moments till the sixth central moments are given below, 

Mi = Ki, 

= K 2 , 

= K 3 

fi 4 = K 4 + 3Kl 

fi 5 = K 5 + 10K 3 K 2 , 

M6 = K 6 + \hK±K 2 + IQKl + 15Jf| . (12) 
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Each term in the expansion for moments in terms of the cumulants is of the form K A * * K^' * such that 
aA + bB + cC + .... = n. 

Also note that, A, B, C... > 1 and 2 < a, b, c, ... < n, where n is the moment that one is interested in. The coefficient 
of any general term in the expansion of the moment in terms of the cumulant is given by, 

n\ 

A\*a\ A * B\*b\ B *C\*c\ c .... ' 
Note that in the figures that appear in the rest of the article we plot the normalized moments defined by, 

II Norm = f^_ (13) 



IV. STATISTICS OF BIPOLAR SPHERICAL HARMONIC COEFFICIENTS 



We classify BipoSH coefficients into four different cases depending upon the form of their characteristic function. 

Case A: h = l 2 ,M = 0, 
Case B: h ^ l 2 ,M = 0, 
Case C: h = l 2 ,M ^ 0, 
Case D: l x ^ l 2 ,M ^ 0. 

To arrive at the distribution of a given BipoSH coefficient, we begin with finding out the characteristic function of 
individual terms involved in the linear combination. Now, the characteristic function of the BipoSH coefficients can 
be written as the product of the characteristic functions of each of the individual terms present in summation. 

This simple scheme works really well for cases A and B. However it only partially works in the the cases C and 
D as there appear terms in the summation which are not statistically independent of each other. For these cases, 
we calculate the moments using the method of characteristic function and then present a general prescription for 
calculating the correction to these moments. All the moments are evaluated in units of standard deviation (a). 



Case A: Bipolar coefficient with ii = I2 = I, M = 



These coefficients are only real, as their imaginary part do not exist. Refer to Appendix [C] for details. In this case, 
all the terms in summation are independent of each other. In the linear combination there will terms with distinct 
distribution functions. Terms with {mi ^ 0,toi = —m 2 } are \ 2 distributed with two degrees of freedom and terms 
with {mi = m 2 — 0} are % 2 distributed with one degree of freedom. For the details of the characteristic function of 
these BipoSH coefficients, refer Appendix IC II 

The n th order cumulant for Aff can be derived to have the following analytical form, 



K n =2"- i (C ; ) n (n-l)!x 



LO \ 



mi 
{mi>0} 



((-i) mi Q™ i; 



Moments for these coefficients can be derived given this form of the cumulants (Eq. [T5|) . We have shown that the 
odd moments for these coefficients oscillate between positive and negative values for even and odd multipoles (I) 
respectively. Example of this behavior can be seen in Fig. [T] and Fig. [2] For coefficients with L ^ 0, mean turns 
out to be zero but rest of the odd moments are non-zero which implies that these coefficients have an asymmetric 
distribution, as seen in Fig. [T] 

A special case of these coeffients are of the form A\f, which are equivalent to the CMB angular power spectrum Cj, 



A 00 
A ll 



(-i) l V2iTTCi 

The characteristic function for these coefficients has the following form, 

-(2i + l)/2 



1 - 



2i(-iycit 
V21TT 



(14) 



(15) 
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FIG. 1: Standard deviation, skewness, kurtosis and 5 th moment of real part of Af®, from 15000 simulations. 
WMAP7 has detected a signal of isotropy violation in these coefficients around the multipole of 200. Hence we 
calculate the statistics of these coefficients up to higher multipoles. 



The Fourier transform of the characteristic function yields the probability distribution function(PDF), which is this 
case is a \ 2 distribution. For even values of multipole (/) the PDF has the form, 



f(x,k) 



(16) 



2 fc/2 a fc/ 1 2 r ( fc / 2 ) a;2 lex P(2a) X - 

0, otherwise 

and for odd values of multipole (I), it is /*(—&, k). Here a = is the related to the isotropic power at multipole I 

and k = 21 + 1 is number of degrees of freedom of the \ 2 distribution. The distribution function for these coefficients 
is asymmetric. The mean for these coefficients are non-vanishing as seen in Fig. [2j this is expected since these are 
the only non- vanishing coefficients under statistical isotropy. 



Case B: Bipolar coefficient with l\ 7^ h, M = 



Even in this case, all the terms in summation are independent of each other. In the linear combination there will 
terms with distinct distribution functions. Terms with {mi 7^ 0,toi = — TO2} are Laplace distributed and terms with 
{mi — m,2 = 0} are modified Bessel of second kind distributed. The details of the characteristic function of real and 
imaginary parts of these BipoSH coefficients can be found in Appendix IC 21 Only even ordered cumulants exist for 
these coefficients, 



k (n=even) ={n-l)\{C h C l2 ) n ' 2 x 



E 



( C LM \' , (r LM s 



mi^0,iri2/0 
mi — — ma 



Note that imaginary part of these coefficients will not have any contribution from the second term in above expression 
for cumulants. Refer Appendix IC 21 for details. Moments of distribution of these coefficients can be obtained given 
the above form for the cumulants(Eq. [T3|) . These coefficients have symmetric PDF, as evident from Fig. |3l 
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FIG. 2: This figure depicts the PDF of some of the coefficients of the form derived from 15000 simulations. 
Coefficients with even multipoles (I) are left skewed and those with odd multipoles (I) are right skewed. 
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FIG. 3: Standard deviation, skewness, kurtosis and 5 th moment of real part of Af® +2 , from 15000 simulations. These 
coefficients have a symmetric PDF. The kurtosis of these coefficients approach that of a Gaussian for high multipoles. 



Case C: Bipolar coefficient with h — I2 — I, M 7^ 



We first calculate the moments of distribution for these coefficients using the characteristic function method assum- 
ing that all terms in the linear combination are independent. In the linear combination for these coefficients there ap- 
pear terms like {mi 7^ 0, 7712 7^ 0} which are Laplace distributed and terms like {mi = 0, ui2 = M},{mi = M, m-i = 0} 
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and {mi — ^2} which are distributed as modified Bessel function of second kind of zeroth order. The details of the 
characteristic function for these coefficients can be found in Appendix IC 31 It is observed that only even ordered 
cumulants exist implying that the distribution of these coefficients is symmetric, 



K( n=eV en) = (tl - 1)!C™ X 



E 

miZm 2 ) (Qmi!m 2 ) ^m 1 m 2 mi 



lm,2 1 



mi>m2 



miVm2-0 



Note that imaginary part of these coefficients will not have any contribution from the last term in above expression 
for cumulants. Refer Appendix IC 31 for details. Moments of distribution of these coefficients can be obtained given 
the above form for the cumulants(Eq. [T3")) . see Fig. @ for illustration. The mismatch in simulations and analytically 
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FIG. 4: Standard deviation and kurtosis of real part of A^f derived from 15000 simulations. The difference between 
corrected and uncorrected analytical moments is prominent at low values of multipole (I). The corrected kurtosis 
can be seen to be in good agreement with the simulation resutls. 



derived moments at low multipoles (I) is due to the assumed underlying independence of the terms contributing to 
the linear combination, which does not hold true for this case. Hence the characteristic function approach does not 
fully characterize the statistics of these coefficients. The moments calculated using the characteristic function method 
need to be supplemented with correction terms, which account for the higher order correlations. Refer to Appendix 
IC 61 for details. However it is found that for certain coefficients the terms involved in the linear combination are all 
independent and the correction term goes to zero, fl are moments calculated using the characteristic function method 
and p, are the corrected moments. 



fi n = /i n + correction . 



(17) 



We show that variance will not have any corrections due the fact that the terms are linearly uncorrclated. However 
kurtosis does have a correction term as seen in Fig. 01 



Case D: Bipolar coefficient with Zi 7^ li,M 7^ 



Similar to the previous case, we begin by finding the moments of distribution for these coefficients using the 
characteristic function method assuming that all terms in the linear combination are independent. In the linear 
combination for these coefficients there appear terms with {mi 7^ 0, ra-i 7^ 0} which are Laplace distributed and terms 
with {mi = 0, VH2 = M}, {mi = M,ni2 = 0} which have modified Bessel function of second kind distribution. The 
details of the characteristic function for these coefficients can be found in Appendix IC 41 Even for these coefficients it 
is found that only even ordered cumulants exist implying that their PDF is symmetric. 



n—even ) 



(n-l)l(C h C l2 ) n / 2 x 



E 



ol— n (siLM 
L l U limi 



mi Vm2 =0 



(18) 



Note that imaginary part of these coefficients will not have any contribution from the last term in above expression 
for cumulants. Refer Appendix IC 41 for details. Just like in the previous case, the moments calculated using the 
characteristic function method are supplemented with correction terms which account for the non-linear correlations, 
see Fig. [5] for illustration. 
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FIG. 5: Standard deviation and kurtosis of real part of AjRf derived from 15000 simulations. The difference 
between corrected and uncorrected analytical moments is prominent at low values of multipole (I). The corrected 
kurtosis can be seen to be in agreement with the simulation results. 



To quantify the agreement between simulations and the analytically derived results we calculate the mean square 
difference. The closeness of fit is seen (Fig. [6| to be inversely proportional to the number of simulations. We observe 
that beyond 10000 simulations good convergence is achieved hence we go up to 15000 simulations to derive all our 
results. 
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FIG. 6: Variation of mean squared difference with number of simulations studied for skewness, kurtosis and 5 th 
moment of the bipolar coefficients. The mean squared difference evaluated for each moment is multiplied with an 
arbitrary constant to bring them to the same scale. The mean squared difference for each moment is found to be 

inversely proportional to the number of simulations. 



Covariance of Bipolar coefficients 

Under statistical isotropy, we show (using Eq. IA6[) that the covariance takes up the following form, 

(At™A%l M ) = C h C h 6 hl{ 6 hl , 2 6 LL , 8 M m> + (-l) h+h+L C h Q 2 5 hV3 S hl{ S LL , <W 

+ ((J h C Vi {-l) h+l '^ '(2h + l)(2Zi + l)S hh 5 VA x 5 lo 6mo6 l , 5 m , } ■ 

We find that the bipolar coefficients are linearly uncorrelated. However this is not a sufficient condition for the 
coefficients to be independent of each other. 

V. DISCUSSION AND CONCLUSIONS 

The temperature field on the CMB sky is usually believed to be a Gaussian random field. Statistical isotropy which 
implies rotational invariance of two point correlation function is an assumption in cosmology and needs to be rigorously 
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tested. Estimators can be constructed aiming at the various kinds of statistical isotropy violations [19( . Knowing the 
PDF of these estimators gives a much better handle on assessing the significance of any statistical isotropy violation 
detection. For example, related analysis has also been carried out to find out the PDF of non-Gaussianity estimators, 
fni, as the significance of a measurement of this parameter depends on knowledge of the full shape of its PDF [20| . 
The two point correlation function is used as a measure of statistics of a Gaussian random field. Correlation function 
is most generally expanded in the BipoSH basis. The coefficients of expansion in this basis encode all the symmetries 
of the correlation function. In this paper we derive the statistical coefficients of these coefficients of expansion. A 
quantitative understanding of the statistics of these coefficients is important, as signal of isotropy violation are being 
searched for in CMB data in these coefficients. 

The strategy has been to calculate the characteristic function for these coefficients and then arrive at the cumulants. 
These cumulants can be easily translated to yield the moments of distribution of the coefficients of expansion. This 
strategy works perfectly well when the terms involved in the expansion of the BipoSH coefficients are independent of 
each other. However we notice that for a certain set of BipoSH coefficients the characteristic function approach works 
only partially. This is due to the presence of non-linearly correlated terms in the expansion of BipoSH coefficients. In 
these cases we give a prescription to account for the contribution of these non-linear correlations to the moments of the 
distribution. In this paper we restrict to calculating the correction to the moments only upto kurtosis, as calculations 
for the correction for higher order moments become increasingly tedious, however the general prescription would work. 

The BipoSH coefficients of the form A®® are directly related to the CMB angular power spectrum. As expected 
these coefficients are shown to have a x 2 distribution with (21 + 1) degrees of freedom using the characteristic function 
method. For rest of the BipoSH coefficients we provide analytical expressions for moments up to any arbitrary order. 
We find that BipoSH coefficients of the form Aff have an asymmetric distribution. The remaining BipoSH coefficients 
are shown to have a symmetric distribution. The BipoSH coefficients of the form Af*j^ (M ^ 0) comprise of terms 
with non-linear correlations amongst them, due to which the analytical moments derived from characteristic function 
method need to be supplemented with a correction term. We give a general prescription to derive the corrections due 
to the presence of these interdependent terms. We explicitly calculate these correction terms only up to kurtosis. All 
these results are tested against extensive simulations. Isotropy violation signal are being sought after in the BipoSH 
representation of CMB maps. A thorough understanding of the statistics of these coefficients is extremely crucial 
to assess the significance any statistical isotropy violation measurement. In the recent past, WMAP7 team claimed 
detection of isotropy violation in V-Band and W-Band maps. This isotropy violating signal was found in the BipoSH 
coefficients Aj® and Aff +2 . The PDF of these coefficients significantly deviate from being Gaussian, particularly, at 
low multipoles. The BipoSH coefficients Aff +2 are found to have a symmetric PDF. 

Interestingly in our study we find that the BipoSH coefficients Af® have an asymmetric PDF, with even multipoles 
(I) being positively skewed and the odd multipoles (I) being negatively skewed. The WMAP team uses band power 
averaged BipoSH coefficients to reduce noise. We find that for full sky and isotropic CMB maps, this averaging results 
in reduced skewness for these coefficients. With experiments like PLANCK it might be possible to achieve similar 
signal to noise ratio for smaller bin sizes. However our study suggests that reducing the bin size, the skewness of these 
coefficients might become considerable. We are currently assessing the implications of these statistics, which is work 
in progress. 
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Appendices 

Appendix A: Statistics of spherical harmonic coefficients 

The temperature fluctuations in the CMB sky maps denoted by AT(h) where h — (6,(f>), can be decomposed in 
the following manner, 

oo m— +1 

AT(h) =J2Y1 a i™Yhn(n) . (Al) 

1=1 m= — l 

where Yi m (n) are the spherical harmonics and ai m are the spherical harmonic coefficients. The expansion coefficients 
can be obtained by taking the inverse tranform of the above equation and can be expressed as, 



a lm = J dn fl Y* n (h)AT(n) . (A2) 
The spherical harmonics can be expressed in terms of the Legendre polynomials, 



Y lm (0,4>) = (~l) m y ^ + m) 7 } 0)6*"* ■ (A3) 

Spherical harmonic coefficients, a; m 's are complex coefficients, 

dim = %lm + iyim, (A4) 

where xi m and yi m are real and imaginary part of a; m and are statistically independent of each other. 

Reality of temperature fluctuations (|Al|) guarantees that the following relation holds for the spherical harmonic 
coefficients, 



-m i 



aim = (-l) m at 

Xlm = (~l) m Xl- 

yim = (-ir + Vm. (A5) 

It is easy to see from the above expressions that when m = 0, the imaginary part of the expansion coefficient 
vanishes. 

CMB temperature fluctuations resulting from the simplest versions of the inflationary paradigm are Gaussian and 
statistically isotropic. The statistical isotropy(SI) takes the form of a diagonal covariance matrix in harmonic space, 

(A6) 

where Ci is the angular power spectrum. For case the of statistical isotropy, the angular power spectrum carries all 
the information about the Gaussian temperature fluctuations. 

The real and imaginary parts of the coefficient a; m , with m ^ 0, are independent Gaussian random variates with 
mean zero and variances given by, 

cr 2 (xi m ) = <J 2 (y lm ) = -Q. (A7) 

However for the coefficients with m = 0, the imaginary part vanishes and the real part are Gaussian random variables 
with mean zero and variance given by, 

a 2 (x lQ )^C,. (A8) 

Appendix B: Characteristic function approach and applications 

The characteristic function of any random variable is defined as the fourier transform of its probability distribution 
function. 

ipx{t) = E[e ltx ] ten. (Bl) 
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Using the characteristic function to arrive at the statistics of random variables forms a very powerful tool, as there 
appear situations in which it is easier to arrive at the characteristic function than the probability distribution function 
itself. This approach is particularly useful in statistical analysis of linear combinations of independent random 
variables. Let Z n be a random variable defined in the following manner, 



(B2) 



where Xi,X%, ...X n are independent random variables which are not necessarily identically distributed and a,'s are 
constants. The characteristic function for such a random variable (Z n ) is given by, 

fz n {t) = <Px 1 ( a i t ) l Px 2 {a2t) tpx n (a„t). (B3) 

Now we discuss a few applications of this technique, which are extensively used in our calculations. 
Let X\ and X 2 be two independent normal variates with zero means and variances a\ and a 2 . The distribu- 
tion of the product of these random variables(Z = X\X 2 ) is given bv[l7l. l2lj. 

fz(z) = (B4) 

7T<7l(72 

where Kq is the zeroth-order modified Bessel function(normal product distribution function). The characteristic 
function corresponding to the above distribution function is given by [22| . 

VzXF, = TT/2 ■ ( B5 ) 



(t 2 + ^ 



Consider the case of linear combination of two normal product distributed random variates. If X\,Y\,X<i,Yi are 
independent Gaussian variates with zero means and variances a\ for X\ , Y\ and a\ for X 2 , Y 2 ■ Then the characteristic 
function of the random variable Z = X\X 2 + Y1Y2 is given by, 

Mt) = IIT W)' (B6) 

where A = l/aia 2 - 

The above characteristic function corresponds to that of a Laplace distribution (Laplace(0, 2A 2 )). 

Another interesting application which is of our interest is that of the difference of squares of two Gaussian 
random variates with zero mean and having the same variance. It is well known that the sum of squares of two 
Gaussian random variates is \ 2 distributed. The difference however is not x 2 distributed, instead it has a modified 
Bessel function of second kind distribution. We demonstrate this using the characteristic function approach, If X 
and Y are two random variables which have normal distribution N(0, a), then X 2 and Y 2 are x 2 distributed and 
their characteristic function is given by, 

^ )= {1-Lh)^ - (B7) 
Using eq. (|B3j) . we obtain the characteristic function for the random variable defined as Z = X 2 — Y 2 , 



(t) m ■ (B8) 

[l + (2(7 2 t) 2 ] 1/2 



Notice that this characteristic function is that of the modified Bessel function of second kind distribution with zero 
order. 



The above illustrated examples are of our particular interest, as they will be used to study the statistics of bipolar 
spherical harmonic coefficients. 



12 



Appendix C: Bipolar statistics 

In order to delve the rich source of information which will be provided by future CMB maps, it is important to device 
methods to detect, isolate and diagnose various possible causes of departure from statistical isotropy. In particular, 
our approach is to look at the statistical behavior of the complex coefficient that arise in bipolar spherical analysis of 
the CMB two point correlation function. 

A^ = <t W +i<™ {1) ■ (CI) 
Owing to the reality of the correlation function, the following relation holds between the BipoSH coefficients, 

a*LM _ r_-i\h+l 2 -L+M aL-M 
hh v I l t l 2 ' 

aLM {r) _ i -,\h+l 2 -L+M aL-MW 
A hh ~ \ L > A hh 

A LM^ = { _ l)( l 1+ l 2 -L + M) + l A L-M^\ (C2) 

The real and imaginary parts of the BipoSH coefficients can be expressed as linear combinations of the elements of 
the covariance matrix in harmonic space, 

aLM^ . V fr, r, - 11, 11, )C LM 
■ n lil 2 ~ / i \- L l\m 1 ^l 2 m 2 <Jl 1 m 1 yi 2 m 2 )^ J l 1 m 1 l 2 m 2 i 

ra\m 2 

A hh = (yi 1 m 1 Xl 2 m 2 + ^(irri! ?/( 2 m 2 )C^m 1 ; 2m2 • (C3) 

ra\ra 2 

The indices in the above expression satisfy the following relations : \l± — l 2 \ < L < li + l 2 and mi + m 2 = M, owing 
to the presence of the Clebsch-Gordon coefficients. 

The BipoSH coefficients can be classified on the basis of the form of their characteristic function, 
Case A: l x = l 2 ,M = 0, 
Case B: h ^ l 2 ,M = 0, 
Case C: h = l 2 ,M ^ 0, 
Case D: l x ^ l 2 ,M ^ 0. 



1. Case A: h = l 2 ,M = 



These BipoSH coefficient can be expressed in terms of the spherical harmonic coeficients of the CMB maps, 



A hh ~ 2—/ a hm 1 Ul 1 m 2 i^i imi i ir 



mim 2 

{(mi^0,m2^0),mi- — m 2 } 



^ a hm i ai 1 m 2 C[ J 1 1 ll i im2 . (C4) 



mi m 2 
{m 1 =m 2 =0} 



We have divided the expansion of the BipoSH coefficients into two type of terms. Terms with {(mi ^ 0, m 2 ^ 0), mi = 
—777,2} and the terms where both mi and 7712 are zero {mi = 777,2 = 0}. This is done because of the fact that each 
of these terms have distinct distributions. The imaginary part of these coefficients vanish owing to the reality of the 
correlation function. 

The real part of these coefficients is given by (Eq. IC3|) , 

A lfh — mi(mi>(l)(^ 1 ) mi ^ :E iimi + yfimi^hmih-mi + ^fiO^OZiO • (C5) 

To arrive at the moments of these BipoSH coefficients, one needs the characteristic function of each term in the 
summation. The first term in the above expression has a x 2 distribution with two degrees of freedom. Its characteristic 
function has the following form (refer Appendix |B|). 
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The second term is x 2 distributed with one degree of freedom and its characteristic function has the following form 
(refer Appendix [Bl . 



<Pz(t) 



Hence, the characteristic function of these BipoSH coefficients is given by, 



n 



i 



„„ L l-(2i(-l)^C^ Il _ mi C Il t)] 

.{mi^O} 



[1 - (2iC$ h0 Cht)) 



1/2 



(C7) 



(C8) 



2. Case B: h =£ l 2 , M = 
The difference in the expansion in this case and the case above is that here l\ ^ fo- 

E r LM , r LM 

u hm 1 ^l 2 m 2 , -^l imi l 2 'm 2 ' / , u hm 1 Ul 2 m 2 { ^l imi l 2 



A 



LM 



(C9) 



mim2 

{(mi/0,m2^0),mi- — 7112} 



mim,2 

{mi— 7712— 0} 



The real and imaginary part of these coefficients are given by the following expressions (Eq. IC3I) , 



11(2 



Hih 



imi x hm2 Vhmi Vl 2 ra 2 )^'l 1 m 1 l 2 m 2 X hO x l 2 oC'l 1 ol 2 Q-> 



m\m 2 

{ (mi ^0,77125^0), mi = — m 2 } 



^ ] {yhm 1 Xl 2 m 2 + x hm i yi 2 m 2 )Ci imi i 2 



(CIO) 



mim 2 

{(mi^0,m2^0),mi--m2} 



Note that the imaginary part of these coefficients does not vanish. 

The first term in the expansion for Af®^ ' and Aj^ ' is Laplace distributed with characteristic function given by 
(refer Appendix IB) . 



<Pz(t) = 



K4 + (C i f mii2m2 i v ^a)2) 



(Gil) 



and the second term in the expansion for Af^f^ has modified Bessel function of second kind distribution with the 
following characteristic function (refer Appendix [Bjl . 



<Pz(t) 



^^ + (ctXi 2m2 tV^ay' 

Hence, the characteristic function for the real part of these BipoSH coefficients is given by, 



(C12) 



PalmM (*) — 



n 



^+Ki lhm2 tv^cr 2 ) 2 ) 

{(m 1 /0,m ! /0),mi=-m 2 } 



n 1 

and the characteristic function for the imaginary part of these BipoSH coefficients is given by, 



(C13) 



¥ a lm(0 (t) 

l l l 2 



n 



^ + C ll2 ™ 2 t v^A) 2 ) 

_{ (mi 7^0,77127^0) ,mi=— ma} 



(C14) 
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3. Case C: h = l 2 ,M ^ 



The expansion of these BipoSH coefficients is split into three parts depending upon the form of the characteristic 
function of each of the terms, 



A LAI 
hh 



E 



On n (~lLM , 



mim2 

{(rai^0,m2/0),mi>m2} 



EOn n r LM J- n n H LM 



m\m 2 
{(mi Vm 2 )=0,mi >m 2 } 



mim 2 
{mi— m 2 } 



The real and imaginary parts of these bipolar coefficients are given by, 

,{R) V 9(t, t, - it, 11, \r LM 

/ j ^\* L l 1 m 1 ^l 1 m 2 yiim 1 yl 1 m 2 )'^l 1 m 1 l 1 i 

mim 2 

{(mi/0,m2/0),rai>m 2 ) 



aLM 1 - 

A hh 



^ t 2xi imi Xl inl2 Ci imi i im2 + ^ ] {xi 1 m 1 x l 1 m 2 yhm i yi 1 m 2 )Ci imi i ir 



mim 2 
{(miVm 2 )=0,mi>m 2 } 



mim 2 
{m 1 =m 2 } 



(C15) 



.4 



LM<'> 
hh 



^ ] 2(yiim 1 x l 1 m 2 + x hin i yi 1 m 2 )Ci imi i im2 + 



mi m 2 

{(mi7^0,m 2 ^0),mi>m 2 } 



E (w 



mim2 
{mi=m 2 } 



(C16) 



The first term in the expansion for -Aj^ ' ano - ^h? 2 ' * s Laplace distributed with characteristic function given by 
(refer Appendix [Bj, 



<PZ® 



(C17) 



The second term in the expansion for Aff ' has a modified Bessel function of second kind distribution. It has the 
following characteristic function (refer Appendix [Bj) , 



<Pz{t) 



(C18) 



The last terms in the expansion for Af ^ and Af® have a modified bessel function of second kind distribution and 
the corresponding characteristic function is given by (refer Appendix [Bj , 



(C19) 



The characteristic function for the real part of BipoSH coefficients can now be easily derived to have the following 
form, 



V a lmW (*) 
ii'i 



n 

*ri\r 

n 



,{(mi^0,m 2 ^0),mi>m 2 } 



mim 
,{(miVm 2 )=0 



' , ,WH(2ca 1B1J at) J 

,mi>m 2 } V 



n i - 



(C20) 
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and the imaginary part of the BipoSH coefficients can be derived to have the following form, 



■(*) 



n 



m\m 2 
{mi/0,m 2 ^0} 
{m 1 >m 2 } 



n 



'2it 



LM 
m\l\m2 



c h tf 



(C21) 



4. Case D: h ^ l 2 ,M ^ 



The BipoSH coefficients in this case will have the following expansion, 



hh ~ / y a £l"»i a Z2«i2 ( -'ZimiZ 2 m 2 + / y a hm 1 U<l 2 m 2 ^l im ih 



mim2 
{mi^0,m 2 ^0} 



mim 2 
j(m 1 Vm 2 )=0} 



The real and imaginary parts of these coefficients can be expressed as fEq lC3[) , 



A hh 



y ] ( x hm 1 Xl 2 m 2 yiim i yi 2 m 2 )C'l 1 m 1 l2 



mim 2 
{mi/0,m 2 /0} 



/ j ■ L llin 1 xi 2 m 2 ^l 1 rnil2m2i 



m\m 2 
j(miVm 2 )=0} 



.4 



hh 



y ] {yiim 1 x l 2 m 2 + x hm i yi 2 m 2 )Ci imi i 2 



mim 2 
{mi^0,m 2 ^0} 



(C22) 



(C23) 



The first term in the expansion for Aj 1 ^ and Af®^ is Laplace distributed and its characteristic function given by 
(refer Appendix iBl) . 



<Pz{t) = 



^(4 + (^iWv^cy 2 ) 



(C24) 



The second term in the expansion for A^^ has a modified Bessel function of second kind distribution. It has the 
following characteristic function (refer Appendix [Bjl . 



(C25) 



Hence, the characteristic function for the real part of BipoSH coefficients is the product of all the characteristic 
functions of individual terms and has the following form, 



11*2 



n 



mi , m2 Ai+(ctx ll2m2 tvchcr 2 ) 2 ) 

,(mi/0,ra 2 /0) 



n * 

M ^z 2 u } H 2+ n A ^v^y 

and the characteristic function for the imaginary part of the BipoSH coefficients is given by, 



l l>2 



n 



,{mi/0,m 2 #0} 



(C26) 



(C27) 
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5. Covariance of BipoSH coefficients 

The unbiased estimator of the BipoSH coefficients is given by, 

^hh ~ ("l)™ 2 aiim 1 «i , 2m2 C^m 1 ; 2 _ m2 • (C28) 

mi 7712 

The covariance of these coefficients is defined in the following manner, 

(M' 1 b 2 A *l^ M ) = \ i~ l ) m2+m2 ^ a l^ a hni2 a \m{ a l'2^0 C ■ (C29) 

Y?7li77l2 TTlim^ J 

The spherical harmonic coefficients (a/ m 's) are Gaussian random variables, hence the four point correlation function 
appearing in the above equation can be written as, 

\ a limi a * 2 m2 a \m' 1 a l' 2 m'J = { a lirn 1 a[ 2m2 ) ( a i';77i' 1 a i2 m 2) 

+ ( a il77lia;*;777' 1 )( a i'2 7772 a '2 m 2) 

+ (ai imi ai' 2m ' 2 )(a^ mfi a* 2m2 ) . (C30) 

Under the assumption of statistical isotropy, the covariance of the BipoSH coefficients can be derived to have the 
following form (Eq |A6p , 

(Aff 2 Afy M ') = C h C lfi (-l) 1 ^ {(2h + l)(2l[ + 1)] 1/2 5 hl2 5 lfil2 S LO S MO S L >oS M ' 

+ C^C^Si^Si^Sll'Smm' + {-^) ll+[2+L Ci 1 Ci 2 8 ll i 2 8 l2Vi 8 LL >5 M M' ■ 

(C31) 

6. Correction to moments due non-linear correlations. 

Consider a random variable define by, 

N 

Z = Y,X t , (C32) 



where X^s are random variable with arbitrary distributions and not necessarily independent and N is total number 
of terms. 

Any arbitrary moment of the distribution of the random variable Z can be expressed as, 



N 



(Z n ) = {(£Xi) n )- (C33) 



In the case where the random variables are all independent of each other, the above expression will acquire this 
simple form, 



N 



(Z")=J2((X l ) n )- (C34) 



However in the case where the random variables are not all independent, the expression for any arbitrary moment 
does not take up the simple form given above. One needs to account for the presence of higher order correlations 
amongst the random variables. This fact needs to be accounted while evaluating each of the moment. 

Specifically while calculating the moments of the BipoSH coefficients we find that the terms appearing in the 
linear combination have non-linear correlations. We evaluate the correction to the moments due to these non-linear 
correlations. We find that there is no correction to the variance as the terms involved in the linear combination turn 
out to be linearly uncorrelated. The corrected kurtosis is derived to have the following form, 

^^4 + (E iY^ )2 . ( C35 ) 
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where second term is the correction term. In the above expression Ki is the cumulant of the i term and Xj and Xj 
are the i th and j th terms in the summation. 

The calculation for correction for higher order moments becomes very tedious, hence we restrict ourselves to 
calculating corrections for moments only upto kurtosis. 
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